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Abstract 


This  contract  funded  a  major  effort  in  bringing  to  maturity,  a  first-principles  theory  of  semiconductor 
lasers  and  amplifiers  based  on  a  fijlly  microscopic  description  of  the  light  semiconductor  material 
coupling.  Significant  breakthroughs  were  made  at  the  microscopic  physics  level  by  explicitly  including 
all  of  the  relevant  excited  bands  in  the  active  region  of  the  device  and  obtaining  quantitative 
agreement  with  gain/index  and  linewidth  enhancement  spectra  measured  for  a  variety  of  Quantum 
Well  structures  at  the  Air  Force  Research  Laboratory,  Albuquerque,  New  Mexico.  The  full  nonlinear 
pde  model  for  high  power  wide  aperture  semiconductor  lasers  was  reformulated  theoretically  so  as 
to  take  advantage  of  speed  and  excellent  parallel  scaling  of  our  ONYX2  in-house  supercomputer. 
Resulting  fi-om  this  is  the  first  interactive  simulator  for  designing  new  generations  of  high  brightness 
semiconductor  devices.  We  were  able  to  begin  building  a  fully  interactive  platform-independent 
optoelectronic  simulator  which  should  have  a  major  impact  on  designing  sophisticated  telecoms 
systems.  Contract  funding  also  allowed  us  to  make  major  inroads  in  explaining  the  novel  creation  of 
an  optically  turbulent  atmospheric  light  guide.  Again  our  excellent  supercomputing  resources 
provided  by  DURIP  awards  allowed  us  to  build  the  sophisticated  parallel  solvers  needed  to  resolve 
extremely  sharp  gradients  associated  with  explosive  intensity  bursts  and  plasma  generation. 
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Bandstructure  and  Microscopic  Many-Body  Gain  Calculations 


Significant  progress  was  made  in  this  task  which  allowed  us  to  include  all  of  the  relevant  valence  and 
conduction  bands  occupied  at  high  carrier  densities  both  within  the  QW(s)  and  the  confining  barrier 
regions.  Inclusion  of  the  barrier  states  was  essential  to  achieve  convergence  in  the  computation  of  the 
linewidth  enhancement  factor  (LWEF),  a  quantity  which  is  the  ratio  of  the  differential  index  to  the 
differential  gain.  This  quantity  is  therefore  an  extremely  sensitive  test  of  the  reliability  of  the 
microscopic  many-body 
theory  unlike  the  gain 
lineshape  which  typically  has 
converged  on  inclusion  of  the 
confined  QW  states  only.  A 
series  of  experiments  carried 
out  by  D.  Bossert  at  the  Air 
Force  Research  Laboratory  in 
Albuquerque,  on  the 
measurement  of  gain,  index 
and  LWEF  spectra  at  difierent 
carrier  densities  provided  a 
testbed  for  the  theory. 

Structures  with  different  QW 
widths  and  different  SCH 
barrier  widths  were  grown  and 
the  above  quantities  measured. 

Results  from  the  theory- 
experiment  comparisons  are  shown  in  the  accompanying  figure.  The  top  picture  is  a  schematic  of  a 
single  InGaAs  QW  Avith  a  narrow  Ga  As  SCH  and  an  AlGaAs  GRINSCH  barrier.  All  of  the 
conduction  and  valence  band  states  included  in  the  combined  bandstructure/many-body  calculation 
are  shown  in  the  top  picture.  The  middle  picture  compares  experimental  (blue  dots)  and  theoretical 
gain  lineshapes  at  different  carrier  densities.  The  bottom  picture  contrasts  the  carrier  density 
dependence  of  the  LWEF  for  a  5  nm  InGaAs  and  a  10  nm  InGaAs  QW.  The  LWEF  for  the  narrow 
well  clamps  with  increasing  carrier  density  whereas  it  shows  a  linear  increase  with  carrier  density  for 
the  wider  QW.  As  high  power  semiconductor  laser  diodes  display  strong  nonuniformity  of  carrier 
density  within  the  device,  we  expect  to  see  a  significantly  different  behavior  for  both  types  of  active 
layer. 


Gain  Medimn  Theory 
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As  an  illustration,  we  ran  a  simulation  of  a  broad  area  device  for  both  QW  structures  at 
levels  and  compared  the  time-averaged  far-field  outputs  for  the  otherwise  identical  lasers, 
field  outputs  are  steering  randomly  on  sub¬ 
nanosecond  time  scales  in  these  devices  due 
to  strong  dynamic  filamentation  instabilities, 
only  time-averaged  quantities  are  meaningful. 

The  figure  on  the  right  shows  that  the  narrow 
QW  structure  (solid  line)  has  an  average  far- 
field  broadening  which  is  about  one  half  of 
the  wider  well  device.  This  result  confirms 
that  our  unique  approach  of  building 
simulation  tools  on  a  predictive  microscopic 
physics  basis  offers  the  possibility  of 
designing  and  optimising  device  performance 
starting  from  the  materials  processing  level. 


high  pump 
As  the  far- 


An  Interactive  Simulator  for  High  Brightness  Diode  Lasers 


Bouyed  by  successes  such  as  the  above,  we  have  embarked  on  an  ambitious  program  to  build  a 
fiilly  interactive  supercomputing  simulator  using  parallel  processing  on  our  Silicon  Graphics 
ONYX2  multi-processor  graphical 
supercomputer.  We  are  designing  the 
interactive  simulator  in  such  a  way 
that  the  user  will  be  able  to  probe  the 
system  and  change  parameters  on  the 
fly.  Pull-down  menus  will  facilitate 
the  utilization  of  tools  to  profile  the 
pump  current,  refractive  index  profile 
and  modify  linear/curved  gratings  for 
wavefront  control  and  surface 
emission. 

A  layout  of  the  simulator  is  displayed 
in  the  picture  to  the  right.  Each  layer 
in  the  structure  indicates  a  particular 
design  feature  or  instantaneous 
snapshots  of  the  relevant  physical 
fields  being  resolved  within  the  active 
layer  of  the  device.  For  example,  it  is 
straightforward  to  add  contacts  of 
general  shape,  grating  sections  etc. 

The  picture  shows  two  contacts 
defining  the  separate  gain  regions  of 
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the  device.  The  narrow  straight  contact  has  gratings  at  both  ends  to  define  a  single  lateral  and 
longitudinal  mode  DBR  master  oscillator.  The  expanding  flare  from  one  end  of  the  master 
oscillator  is  the  power  amplifier  section.  The  facet  reflectivities  at  the  wide  end  of  the  flare  can  be 
varied  and  one  usually  assumes  an  AR  coating  with  reflectivities  as  low  as  0.0005.  Even  with  such 
low  reflectivity,  the  device  can  readily  switch-on  as  a  multi-mode  laser  when  the  master  oscillator 
current  is  modulated.  The  layer  beneath  the  contacts  shows  a  contour  of  the  instantaneous  carrier 
density  distribution  thoughout  the  active  layer  region.  The  dark  troughs  at  the  edges  of  the 
expanding  flare  represent  large  unsaturated  carrier  densities  and  these  are  a  major  source  of 
degradation  of  the  far-field  of  this  device.  By  adding  current  profiling  and  cavity  spoilers,  such 
devices  can  be  pushed  to  higher  output  powers  while  retaining  high  brightness.  This  simulation 
tool  is  becoming  an  invaluable  resource  for  future  high  brightness  system  design  and  optimization. 
A  specific  illustration  of  a  high  brightness  laser  design  which  utilizes  current  profiling  and  a  cavity 
spoiler  is  shown  in  the  picture  on  the 
right.  The  top  left  frame  shows  a  the 
carrier  density  profile  in  a  3-contact 
bow-tie  laser  with  a  flat  current  pump 
profile  throughout.  The  narrow 
central  section  is  meant  to  act  as  a 
spatial  mode  filter  and  has  a  grating 
section  for  longitudinal  mode 
filtering.  A  current  of  100  mA  is 
applied  to  this  narrow  contact  and 
each  flare  has  a  current  of  1.5  A.  The 
lower  picture  on  the  left  shows  the 
identical  device  but  now  with  a 
modified  current  profile  and  a  cavity 
spoiler  placed  at  the  center  of  the  narrow  section.  The  shapes  of  the  current  profiles  are  shown  to 
the  right  of  each  picture.  The  dark  stripes  at  the  edge  of  the  flares  in  the  top  picture  indicate  the 
formation  of  regions  of  high  unsaturated  carrier  density.  These  promote  filamentation  instabilities 
in  the  device  and  lead  to  significant  far-field  degradation.  In  addition  the  fuzzy  areas  outside  the 
straight  contact  edges  are  regions  where  the  optical  intensity  has  bleached  out  the  carriers  and  the 
gain  material  is  transparent  in  these  regions.  Increasing  the  pump  current  leads  to  stronger 
bleaching  and  an  increased  tendency  for  filamentation.  The  profiled  contact  suppresses  the  growth 
of  unsaturated  carriers  at  the  contact  edges  and  the  cavity  spoiler  acts  as  a  lossy  element  for  the 
bleached  regions  of  the  device.  Consequently,  the  brightness  of  the  lower  device  is  considerably 
enhanced.  This  build-up  of  unsaturated  carriers  has  been  seen  experimentally  by  the  group  of  J. 
Mclnemey  (Cork)  in  flared  contact  devices  by  imaging  the  spontaneous  emission  through  a 
transparent  contact  on  the  surface  of  the  device.  This  group  has  used  an  older  version  of  our 
device  simulation  software  to  design  profiled  contacts  and  have  demonstrated  improved  far-field 
emission  in  broad  area  devices  [  P.O’Brien,  P.M.W.  Skovgaard,  J.G.  Mclnemey  and  J.S.  Roberts, 
Electron.  Letts.,  34,  1943  (1998);  P.M.  W.  Skovgaard,  P.  O'Brien  and  J.G.  Mclnemey,  Electron. 
Letts.,  34,  1950  (1998)] 
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An  Optoelectronic  Systems  Simulator 


The  gain  spectra,  once  calculated  for  a  specific  structure,  can  be  used  as  the  active  layer  in  a  host 
of  different  geometry  laser  structures  including  high  power  edge/surface  emitters,  low  power 
VCSEL:s  and  low  power  edge  emitters  typically  employed  in  telecommunications  systems.  In 
telecommunications  systems,  the  semiconductor  lasers  have  strong  lateral  and  transverse  optical 
confinement  and  can  be  modelled  as  1+1  D  systems.  This  allows  us  to  quantitatively  resolve  the 
multi-longitudinal  mode  behavior  of  a  general  device  in  a  computationally  fast  and  efficient 
manner.  The  semiconductor  laser  solver  can  now  become  an  individual  module  in  a  sophisticated 
optical  systems  simulator  and  various  other  optical  components  in  realistic  telecommunications 
systems  ( DFB/DBR  lasers,  fibers,  AR  and  HR  coatings,  beam  splitters,  external  cavities,  etc)  can 
be  built  and  accessed  via  pull-down  menus  within  an  optical  systems  graphical  GUI  editor.  We 
embarked  building  this  simulator  at  the  beginning  of  the  Fall  semester  1999  and  are  making  rapid 
progress  in  assembling  the  simulator.  As  the  simulator  is  again  built  on  a  fijlly  microscopic  basis,  it 
should  become  a  powerful  interactive  simulation  tool  for  designing  conventional,  soliton-based  or 
novel  chaos-based  telecommunications  systems.  Our  plan  is  to  build  a  platform  independent 
simulation  tool  which  will  take  advantage  of  the  rapid  increase  in  CPU  clock  speeds  of  new 
generations  of  intermediate  end  graphical  workstations.  Ideally  the  simulator  will  require  at  least  a 
dual  CPU  box  with  a  high  bandwidth  communications  interface  so  that  the  heavy  parallel 
graphical  diagnostics  (intensity  traces  without  and  with  detector  averaging,  power  spectra,  eye- 
diagrams  etc)  can  be  off-loaded  from  the  computationally  heavy  core  simulator.  A  rather 
primitive  but  working  version  of  a  simulator  window  is  shown  below. 


Fast  Ektemal  Ftekf  Slow  BdemRl  Field 


~r“ .  '  I" . . . . . — - - 

299^0  W9J;  300.0 

ftwjpTHiJ 


6 


The  top  frame  shows  the  z-resolved  (laser  axis)  forward  (red),  backward  (green)  internal  field 
intensities  and  the  carrier  density  (blue)  along  the  device.  The  left  middle  frame  shows  the  finely 
resolved  evolution  of  the  laser  output  and  the  right  middle  frame  a  detector  averaged  response 
over  a  much  longer  time  window.  The  bottom  frame  depicts  the  time  evolution  of  the  spectrum  of 
the  output  field.  The  case  depicted  here  is  a  short  (single  mode)  laser  cavity  with  an  external 
reflector  with  a  reflectivity  of  1%  placed  about  1  m  away  from  the  original  cavity.  The  feedback 
induces  an  instability  in  the  device  which  manifests  itself  as  a  chaotic  pulsing  of  the  output  (left 
middle  frame)  with  pulse  durations  of  about  15-20  ps.  The  coarse,  detector  averaged  output 
shows  the  characteristic  power  dropouts  of  the  device  and  the  spectrum  the  extremely  broad 
chaotic  output  spectrum. 


High  Power  Femtosecond  Light  Strings 


Our  contract  has  enabled  us  to  leverage  a  preliminary  study  of  a  very  promising  new  field  which 
has  the  potential  of  providing  a  broad  base  of  important  military  and  civilian  applications 
opportunities.  The  ACMS 
group  were  the  first  to 
identify  the  highly  dynamic 
mechanism  by  which  a  high 
power  femtosecond- 
duration  laser  pulse  can 
form  a  nonlinear  waveguide 
in  air  [98-2].  The  critical 
power  for  self-focusing  in 
air  is  approximately  2  GW 
and  an  ultrashort  laser  pulse 
with  incident  powers 
exceeding  this  level  can 
critically  focus  in  air, 
causing  optical  breakdown 
at  where  the  intensity 
becomes  extremely  high  (  « 

10^^  -  10*'*  W/cm^).  A  narrow  plasma  filament  is  created  by  the  intense  pulse  with  diameter 
ranging  between  SOpm  (IR)  to  120iam  (UV)  and  of  the  order  of  a  meter  or  longer  in  length.  We 
identified  the  underlying  physical  mechanism  for  the  formation  of  this  light  guide  as  a  form  of 
“dynamic  spatial  replenishment”  whereby  the  pulse  can  undergo  critical  collapse  (self-focusing) 
more  than  once  within  a  single  pulse.  This  dynamic  mechanism  is  illustrated  in  the  four  frames  in 
the  picture  on  the  right.  The  computation  here  was  carried  out  in  radial  symmetry  and  the  radial 
coordinate  corresponds  to  the  transverse  coordinate  in  each  frame.  The  coordinate  along  the  pulse 
is  the  time  coordinate.  The  pulse  (top  picture)  first  begins  to  self-focus  on  the  leading  edge  and 
the  collapse  singularity  is  regularized  by  the  generation  of  a  narrow  plasma  filament  which  is  left 
in  the  wake  of  the  pulse  i.e  the  plasma  recombines  on  a  ps  time  scale.  The  narrow  plasma  filament 
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acts  as  a  strongly  defocusing  lens  which  counteracts  the  self-focusing  and  strongly  evacuates  the 
light  immediately  behind  the  initial  focus.  The  plasma  turns  off  as  the  intensity  drops  and  a 
secondary  collapse  can  occur  if  the  power  behind  the  leading  focus  exceeds  critical.  This  process 
can  recur  a  number  of  times.  Each  occurrence  leaves  a  narrow  plasma  filament  in  the  wake  of  the 
intensity  focus  and  these  filaments  recombine  on  a  much  slower  scale. 

When  the  power  is  cranked  up  such  that  the  incident  laser  pulse  contains  tens  to  hundreds  of 
critical  powers  the  situation  becomes  extremely  interesting.  For  example,  an  experiment  by  a 
group  in  Jena,  Germany  using  a  2.5  TW  100  fs  laser  pulse  were  able  to  probe  the  upper 
atmosphere  (  « 12  km)  by  delivering  a  white  light  continuum  spectral  probe  from  which  they  were 
able  to  detect 
absorption  lines  fi'om 
various  molecular 
species.  In  this 
experiment,  the 
incident  pulse 
contains  on  the  order 
of  1000  critical 
powers.  The  key 

physics  associated  with  this  process  is  that  femtosecond-duration  pulses  cause  optical  breakdown 
primarily  via  an  instantaneous  multi-photon  rather  than  a  cumulative  avalanche  process.  Funding 
in  the  present  contract  helped  us  develop  a  sophisticated  adaptive  mesh,  parallel  algorithm  which 
enabled  us  to  drop  the  radial  symmetry  assumption  above  and  extend  our  study  to  3D  +  time  [99- 
15].  Although  we  could  not  hope  to  simulate  the  Jena  experiment,  we  were  able  to  take  pulses 
with  tens  of  critical  powers  and  demonstrate  a  fundamentally  new  physical  phenomenon  -  the 
creation  of  optically  turbulent  light  strings.  The  accompanying  picture  shows  three  separate  time 
slices  of  the  optical  pulse  and  displays  the  chaotic  creation  of  intensity  “hot-spots”  both  across  ( 
via  a  modulational  instability)  and  along  (  generalization  of  the  radially  symmetric  “dynamic 
spatial  replenishment”  scenario)  the  pulse.  Accompanying  each  light  intensity  focus  is  a  plasma 
channel. 


Air  Force  Links 


Important  collaborative  links  to  the  Air  Force  Research  Laboratory  in  Albuquerque  enabled  us  to 
fully  test  the  microscopic  many-body  theory  against  structures  that  were  completely  characterized 
at  the  growth  level.  It  is  important  to  stress  that  this  link  with  D.  Bossert  was  crucial  to  our 
success  in  anchoring  the  combined  bandstructure/many-body  calculations  as  the  growth  details  of 
most  materials  are  not  available  to  us  and  considered  proprietary  by  industry.  The  experimentally 
measured  gain,  index  and  linewidth  enhancement  factor  spectra  on  these  QW  devices  enabled  us 
to  identify  the  mechanism  behind  the  LWEF  clamping. 

Our  group  has  also  established  direct  links  to  the  DELO  branch  ( C.  Denman  and  A.  Gavrielides) 
on  two  projects.  We  have  had  preliminary  meetings  on  modeling  double  clad  optical  fibers  for 


8 


high  brightness  applications  and  this  area  has  now  been  identified  as  a  priority  under  our  new 
contract  award.  Studies  on  chaotic  dynamics  of  semiconductor  lasers  with  external  cavities  were 
also  undertaken  by  us  [99-2]  and  we  identified  the  important  role  of  isolated  laser  longitudinal 
modes  in  offering  a  WDM  chaotic  synchronization  and  communications  paradigm  for  a 
telecommunications  system.  We  are  currently  in  communication  with  A.  Gavrielides  at  DELO  and 
plan  to  use  our  optoelectronic  simulator  discussed  above  to  implement  a  full  chaotic 
communications  interactive  testbed. 


Industrial  Collaborations  and  Links 

The  present  contract  leveraged  initial  contacts  with  Opto  Power  Corporation  in  Tucson  and  led  to 
the  award  of  an  NSF  GOALI  grant  to  work  directly  with  this  company.  We  are  currently  pursuing 
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Abstract 

We  report  on  the  design  and  fabrication  of  tapered  cavity  grating  coupled  surface- 
emitting  DBR  lasers  in  the  980-nm  regime.  A  curved  second  order  grating  is  used  at  the 
end  of  a  tapered  gain  section  to  provide  feedback  as  well  as  collimated  surface  out- 
coupling.  A  detailed  numerical  analysis  shows  a  far-field  central  lobe  power  of  75%  is 
possible  up  to  twice  threshold.  Near  diffraction-limited  collimated  surface-emitting 
output  with  moderate  power  of  about  150  mW  is  obtained  under  continuous  operation. 
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I. 


INTRODUCTION 


High  brightness  semiconductor  lasers  are  important  light  sources  for  applications  such  as 
pumping  of  EDFAs,  Raman  amplifiers,  and  for  free  space  optical  communications.  In 
order  to  obtain  high  power  from  a  semiconductor  laser,  a  large  aperture  is  needed  to 
avoid  the  problems  of  optical  damage  and  high  current  density  operation.  However,  wide 
aperture  lasers  undergo  rapid  filanientation  and  multi-lateral  mode  operation.  To  prevent 
multi-lateral  mode  oscillation,  some  method  of  mode  discrimination  is  needed.  To  date, 
the  most  reliable  concepts  have  been  the  master  oscillator  power  amplifiers  (MOPAs)  and 
unstable  resonator  lasers.  High  power  single  mode  operation  has  been  reported  in  both  of 
these  edge-emitting  devices  [l]-[3].  However,  due  to  their  high  degree  of  anamorphicity 
and  astigmatism  in  the  output  beam,  coupling  of  these  lasers  to  a  single-mode  waveguide 
or  fiber  is  problematic.  Highly  specialized,  high  numerical  aperture  optics  are  required 
for  efficient  coupling  into  single-mode  fibers. 

To  alleviate  the  need  for  external  optics  and  related  optomechanical  alignment 
problem,  grating-coupled  surface-emitting  lasers  (GCSEL)  can  be  used.  In  comparison  to 
edge-emitting  lasers,  surface-emitting  devices  have  several  advantages:  low  divergence 
angle,  compatibility  with  beam  forming  elements,  single-longitudinal  mode  operation,  on 
wafer  testing  and  2-D  array  integration.  Based  on  the  feedback  grating  configurations, 
GCSELs  can  be  divided  into  two  types:  the  first  uses  first-order  grating  for  feedback  and 
detuned  second-order  grating  for  out-coupling,  the  second  uses  a  single  second-order 
Bragg  grating  for  both  feedback  and  out-coupling.  The  first  configuration  separates  the 
out-coupler  from  the  oscillator  and  gives  the  maximum  flexibility  [4].  However,  because 
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feedback  and  out-coupling  have  different  grating  periods,  they  may  experience 
wavelength  detuning  and  beam  steering  as  operating  condition  changes.  The  second 
configuration  combines  the  feedback  and  out-coupling  into  a  single  second-order  Bragg 
grating.  While  the  beam-forming  capability  is  somehow  limited  in  this  configuration, 
stable  normal  to  the  surface  out-coupling  is  intrinsically  guaranteed  without  any  beam 
steering  [5].  In  addition,  by  matching  the  curvature  of  the  second  order  grating  to  the 
phase  front  of  the  forward  propagating  wave,  collimated  outcoupling  beam  can  be 
achieved.  Finally  due  to  the  fact  that  the  emitting  aperture  is  significantly  larger  than  the 
edge-emitting  laser,  catastrophic  optical  damage  can  also  be  greatly  reduced. 

Here  we  study  the  design  of  a  tapered  cavity  surface-emitting  laser  with  a  curved 
grating.  The  paper  is  organized  as  follows:  in  section  II  we  present  the  numerical  analysis 
based  on  a  comprehensive  semiconductor  laser  modeling.  Section  III  describes  the  design 
and  fabrication  of  the  laser.  In  section  IV  the  experimental  results  are  described.  We  then 
present  discussion  and  concluding  remarks  in  section  V. 

II.  MODELING 

The  schematic  illustration  of  the  tapered  cavity  surface-emitting  DBR  laser  is  shown  in 
Fig.  1 .  The  laser  consists  of  an  index-guided  ridge  waveguide  followed  by  a  gain-guided 
tapered  section.  A  curved  second-order  Bragg  grating  provides  feedback  and  surface- 
emission  while  the  other  facet  can  be  defined  by  cleaving.  The  results  obtained  with  the 
tapered  DBR  resonator  shown  in  Fig.  1  will  depend  strongly  on  the  curvature  and  in¬ 
plane  reflectivity  of  the  second-order  grating  used  on  the  flared  end  of  the  resonator.  To 
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obtain  a  narrow  far-field  pattern,  for  example,  this  grating  must  be  designed  to  match  the 
wavefront  curvature  of  the  beam  from  the  ridge  waveguide  on  the  far  end  of  the  flare, 
thereby  requiring  a  semi-stable  resonator  configuration.  Unfortunately,  the  strength  of  the 
resulting  counter-propagating  field  increases  the  tendencies  toward  fi  lamentation 
typically  found  in  broad-area  semiconductor  devices  [6]-[9]. 

These  tendencies  can  be  strongly  attenuated  with  an  unstable-resonator  design. 
First  proposed  for  Fabry-Perot  cavity  [10]  [11],  these  devices  achieved  brightness  levels 
on  the  order  of  20-30  times  greater  than  that  for  stable  cavities.  Typical  magnifications 
for  such  performance  are  approximately  2  to  3  [12],  with  greater  magnifications  leading 
to  larger  losses  and  therefore  lower  efficiencies.  O’Brien  extended  these  concepts  to  the 
case  of  an  edge-emitting  flared  unstable  resonator  using  straight-line  distributed-Bragg- 
reflectors  (DBRs)  to  define  the  cavity  [13].  As  shown  below,  the  mismatch  between  the 
grating  and  field  curvatures  results  in  extremely  large  far-field  angles  if  this  unstable 
resonator  is  used  in  a  surface-emitting  geometry.  Thus  our  design  is  a  compromise 
between  a  filament-resistant  unstable  cavity  and  the  potential  for  a  diffraction-limited 
beam  for  lower-power  operation  with  the  semi-confocal  curved-grating  cavity. 

In  this  section,  we  present  numerical  results  showing  the  sensitivity  of  far-field 
beam  size  to  power  level.  Our  approach  is  to  first  find  a  grating  curvature  that  produces  a 
narrow  far  field  at  a  specific  pumping  current;  we  then  increase  the  power  level  with  the 
expectation  that  carrier-induced  refractive  index  changes  will  modify  the  power  in  the 
side  lobes.  Our  model  has  been  described  in  detail  previously,  and  includes  the  effects  of 
counter-propagating  fields,  lateral  diffraction,  carrier  diffusion,  and  a  comprehensive  gain 
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model  that  incorporates  many-body  effects  [14]  [15].  Gain/index  coupling  via  the 
linewidth  enhancement  factor  a  is  inherent  in  the  gain  model;  no  explicit  input  of  a 
spatially  averaged  factor  is  required.  Including  curved-grating  effects  and  neglecting  gain 
and  absorption  features  described  previously  [15],  the  forward  and  backward-propagating 
field  envelopes  F{x,z,t)  and  B{x,z,i)  in  the  grating  region  are  given  by 

J_  ^  dF^  ^  (1 ) 

V,.  dz  2k  dx^  \  •  I  \ 


ISB^ _  aB{^ ^  j_d^B(x  z,,) _ (2) 
dz  2k  dx^  \  ’  j  \  !  \  2  , 


v^,  dl 


where  x  and  z  are  the  lateral  and  longitudinal  coordinates,  v  is  the  group  velocity. 


k  =  con  /c  is  the  magnitude  of  the  field  wavenumber,  k  is  the  in-plane  coupling 

coefficient,  e  is  the  outcoupling  coefficient  for  surface  emission,  and  7?(z)  is  the  z  - 
dependent  radius  of  curvature  for  a  concentric  circular  grating.  Note  that  the  quadratic 
phase  factor  applies  only  in  the  paraxial  regime  typical  of  the  ten  degrees  flare  angle 
found  in  most  tapered  devices.  Also  note  that  the  carrier  density  and  electronic 
polarization  equations  are  not  affected  by  the  curved  grating. 

To  obtain  far-field  intensity  profiles,  we  first  propagate  Eqs.  (1)  and  (2)  for  about 
10  nanoseconds,  allowing  time  for  oscillations  to  settle.  We  then  obtain  a  near-field 
amplitude  profile  Efj{x,z^  consisting  of 


Ef,  (x,  z)  =  F{x,  zy'*'^' -H  B(x, z)e 


ikx-‘/2K(z) 
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where  the  phase  of  F{x,z)  and  D{x,z)  are  corrected  for  the  way  in  which  the  grating 
out-diffracts  the  in-plane  wavefront.  The  far-field  intensity  is  then  obtained  from  the 
magnitude-squared  of  the  Fourier  transform  of  {x,  z) . 

The  far-field  intensity  outcoupled  from  a  linear  grating  at  the  end  of  the  flare  is 
shown  in  Fig.  2.  For  this  geometry,  R{z)  =  co  for  all  z ,  so  no  correction  occurs  in  Eqs. 

(1)  and  (2).  The  resulting  phase  mismatch  between  the  field  and  grating  curvature  results 
in  an  enormous  far  field  in  comparison  with  the  0.38  degree  full-width  at  half-maximum 
(FWFIM)  diffraction  limit  for  a  150  pm  aperture.  The  simulation  assumes  k  -  e  -\0cm~' 
for  the  coupling  coefficients;  the  gain  medium  is  a  10  nm  InGaAs  quantum  well  (QW) 
sandwiched  between  graded  index  AlxGai-xAs  layers,  where  x  varies  from  0.1  to  0.6  over 
a  distance  of  85  nm. 

For  a  circular  curved  grating,  on  the  other  hand,  with  an  inner  radius  that 
approximately  matches  the  length  of  the  flare.  Fig.  3  shows  that  a  near-diffraction-limited 
spot  size  with  reasonable  side  lobes  is  obtained  close  to  threshold  (1.03  Ju,)-  At  higher 
power,  we  expect  the  a-factor  to  change  the  wavefront  curvature  from  its  near-threshold 
value,  thus  affecting  the  far-field  profile.  Simulations  show  that,  excluding  thermal 
effects,  this  does  not  become  important  until  about  twice  threshold.  Above  this,  the 
temporal  oscillations  do  not  settle.  Dynamic  instabilities,  possibly  of  the  type  described 
by  Lowery  [16],  are  observed  numerically.  As  shown  in  Fig.  4,  the  on-axis  beam  angle  at 
a  particular  time  can  retain  it's  low-divergence  character;  unfortunately,  the  increase  in 
side-lobe  intensity  decreases  the  fraction  of  power  in  the  central  lobe  (or  power  in  the 
bucket)  to  67%  (from  80%  in  Fig.  3).  The  asymmetry  in  the  far  field  is  indieative  of  the 
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temporal  variations  at  this  power  level  (2.05  Jti,)-  Note  that  at  even  higher  powers,  the 
resulting  far-field  intensity  is  significantly  degraded  (Fig.  5)  and  again  varying  with  time. 

To  summarize:  a  curved  grating  at  the  end  of  the  flare  is  essential  for  a 
diffraction-limited  spot.  Assuming  a  reduction  in  central  lobe  power  to  about  75%  and  a 
small  DFB  outcoupling  coefficient  is  acceptable  for  a  given  application,  a  reasonable 
power  variation  is  then  possible.  Given  the  huge  design  “space”,  including  pumping 
current,  radius  of  the  curved  grating,  grating  shape  (which  we  have  not  considered  here), 
outcoupling  and  in-plane  DFB  coupling  coefficients,  and  linewidth  enhancement  factor, 
our  simulations  show  the  possibility  of  obtaining  a  narrow-angle,  surface-emitting  beam 
up  to  about  twice  threshold. 

III.  DESIGN  AND  FABRICATION 

The  3 -pm  wide  ridge  is  designed  to  only  support  the  fundamental  lateral  mode.  The  ridge 
length  is  400-pm.  Close  to  the  junction  between  the  ridge  and  the  tapered  section  are  two 
spoiling  grooves  located  at  the  two  sides  of  the  ridge.  The  function  of  the  spoiling  groove 
is  to  prevent  any  feedback  from  the  tapered  section  to  the  unpumped  region  in  the  ridge 
side.  Therefore  the  ridge  and  spoiling  groove  combination  functions  like  a  spatial  mode 
filter,  which  will  promote  the  fundamental  mode  and  suppress  high  order  lateral  modes. 
Follow  the  ridge  is  the  tapered  section,  which  provides  most  of  the  gain  needed  for 
lasing.  The  length  of  the  tapered  section  is  1 800-pm,  and  the  end  width  of  the  taper  is 
225-pm.  The  taper  has  a  full  angle  of  about  7  degree,  which  is  designed  to  accommodate 
the  free  diffraction  of  the  beam  inside  the  tapered  section.  The  resonator  consists  of  a 
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cleaved  mirror  on  the  ridge  side  and  a  curved  second-order  Bragg  grating  on  the  tapered 
side.  The  cleaved  mirror  is  HR-coated  to  decrease  the  lasing  threshold  and  to  increase  the 
out-coupling  efficiency  on  the  grating  side.  The  second-order  grating  provides  both 
feedback  and  surface  out-coupling.  The  curvature  of  the  grating  is  designed  to 
approximately  match  the  phase  front  of  the  forward  propagating  beam,  in  this  case,  the 
center  of  curvature  is  located  at  1800  pm  from  the  tapered  end.  This  will  collimate  the 
output  beam  so  that  it  has  a  small  divergence  angle. 

The  devices  were  fabricated  from  a  strained  InGaAs/GaAs  three-quantum  well 
GRINSCH  structure.  The  layer  structure  consists  of  a  1.5-pm  thick  n-AlxGauxAs 
(x~60%)  cladding  layer,  a  0.12-pra  thick  undoped  AlxGai.xAs  (x~60%-25%)  GRIN 
optical  confinement  layer,  three  6-nm  InGaAs  quantum  wells  with  20-nm  GaAs  barrier,  a 
0.12-pm  thick  AlxGai.xAs  (x~60%-25%)  GRIN  optical  confinement  layer,  a  0.3-pm  thick 
p-AlxGai.xAs  (x~60%)  cladding  layer,  a  10-nm  GaAs  etch-stop  layer,  0.9-pm  thick  p- 
AlxGai-xAs  (x~60%)  cladding  layer,  and  a  0.1-pm  thick  p'^-GaAs  cap  layer. 

First,  the  top  contact  (Ti/Pt/Au)  for  the  ridge  and  tapered  section  is  formed  by  lift¬ 
off  process.  The  ridge  section  and  grating  section  are  simultaneously  etched  down  by 
electron-cyclotron-resonance-reactive-ion  etching  (ECR-RIE).  The  etch-depth  is 
controlled  by  timing  to  ensure  that  the  etching  stops  about  0.1 -0.2  pm-  above  the  etch- 
stop  layer.  A  selective  wet  etching  is  then  used  to  remove  the  remaining  AlGaAs  down  to 
the  GaAs  etch-stop  layer.  Here  the  inclusion  of  an  etch-stop  layer  is  important.  It  allows 
us  to  perfectly  control  the  fabrication  of  the  ridge  so  the  lateral  mode  profile  can  be 
ensured.  It  also  enables  us  to  control  the  position  of  the  grating  with  great  accuracy.  After 
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the  grating  area  is  etched  down  to  the  etch-stop,  curved  second  order  Bragg  gratings  are 
defined  by  electron  beam  lithography.  Since  the  grating  section  is  identical  to  the  gain 
section  (no  regrowth),  the  grating  period  is  designed  so  the  lasing  wavelength  is  slightly 
red-shifted  from  the  gain  peak.  This  is  necessary  to  reduce  absorption  losses  in  the 
passive  section.  This  choice  however  makes  the  lasing  less  efficient.  It  also  increases  the 
a-factor,  which  may  cause  a  tendency  towards  filanientation.  The  gratings  are  then 
etched  by  a  second  ECR-RIE.  The  etch-depth  of  the  grating  is  about  150  nm.  After  that 
spoiling  grooves  are  defined  by  wet  etching.  The  etching  is  timed  to  etch  through  the 
active  region  to  prevent  any  unwanted  feedback.  The  device  is  then  lapped  in  the 
backside  down  to  about  100-pm  total  thickness.  A  Ni/Ge/Au  n-side  metallization  is  then 
deposited.  The  device  is  annealed  using  rapid  thermal  annealing  (RTA)  to  form  an  ohmic 
contact.  The  device  is  cleaved  on  the  ridge  side  and  the  mirror  is  coated  for  a  reflectivity 
of  about  90%.  Finally,  the  device  is  bonded  p-side  up  on  a  gold  plated  14-pin  mount  and 
wire-bonded  for  testing. 


IV.  EXPERIMENTAL  RESULTS 

The  device  is  tested  at  room  temperature  under  pulsed  and  continuous  operation.  The 
laser  is  placed  p-side  up  on  a  copper  plate  whose  temperature  is  controlled.  Since  the 
laser  is  die-bonded  p-side  up,  the  thermal  resistance  is  expected  to  be  very  high.  For  the 
pulsed  operation,  the  pulse-width  is  fixed  at  5-ps,  and  the  duty-cycle  is  changed  to 
investigate  the  thermal  effect.  The  L-I  characteristic  under  pulse  operation  with  duty- 
cycle  of  1%  and  10%  are  shown  in  Fig.  2.  There  is  little  difference  between  these  two 
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duty-cycles.  Under  the  pulsed  operation,  the  threshold  current  is  650  niA  (corresponding 
to  a  threshold  current  density  of  315A/cni^)  and  the  external  differential  quantum 
efficiency  is  24.4%.  Output  power  of  120  mW  is  obtained  under  the  pulsed  operation  at  a 
pump  current  of  1  A,  which  is  the  limit  of  our  pulse  current  source.  There  is  no  sign  of 
degradation,  which  indicates  that  higher  power  is  possible  at  higher  current.  The  laser  is 
then  tested  under  continuous  operation  with  p-side  up.  The  L-I  curve  under  CW  is  also 
shown  in  Fig.  2.  The  threshold  current  is  700  niA  (corresponding  to  a  threshold  current 
density  of  339A/cm^)  and  the  external  differential  quantum  efficiency  is  15.5%.  Output 
power  of  150  mW  is  obtained  under  CW  operation  at  a  pump  current  of  1.5  A.  The  peak 
output  power  is  limited  by  the  thermal  rollover  as  the  current  goes  higher.  The  drop  of  the 
efficiency  and  power  under  CW  compared  to  pulsed  operation  is  also  observed  in  other  p- 
side  up  mounted  broad-area  lasers.  This  is  mostly  due  to  the  large  thermal  effect  under 
CW  operation. 

To  investigate  the  beam-quality  of  the  device,  its  far  field  at  different  driving 
currents  is  measured.  The  far  field  is  measured  by  directly  recording  the  beam  on  a  CCD 
camera  200  mm  away  from  the  near  field.  This  measurement  setup  does  not  correct  any 
residual  phase  curvature  that  may  be  present  in  the  near  field,  so  this  method  measures 
the  "true"  divergence-angle  of  the  far  field.  The  measured  far-field  patterns  under  CW 
operation  at  different  currents  are  shown  in  Fig.  3.  Below  1.3  A  (1.9  lu,),  the  far  field  had 
only  one  spot.  As  the  current  increases  from  the  threshold  up  to  1.1  A  (1.6  In,)  the 
divergence  angle  of  the  beam  stays  below  0.25”(FWHM),  which  corresponds  to  1 .1  times 
diffraction  limited  value  for  a  225  pm  wide  grating.  As  the  current  increases  further,  the 
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divergence  angle  increases  slowly.  The  divergence  angle  increases  to  0.30°  at  1.3  A, 
which  coiresponding  to  approximately  1.4  times  diffraction  limited.  As  the  current 
inci eases  beyond  1 .3  A,  a  second  lobe  starts  to  appear  and  the  far  field  does  not  maintain 
a  single  lobe.  We  believe  the  increase  of  the  divergence  angle  is  mostly  due  to  the 
thermal  lensing  effect.  More  efficient  heat  sinking  and  better  packaging  should  increase 
the  range  of  the  current  for  near  diffraction-limited  operation. 

V.  CONCLUSION 

We  have  carried  out  a  detailed  numerical  analysis  of  a  surface-emitting  laser  with  curved 
gratings.  The  simulation  results  show  that  a  near-diffraction  limited  collimated  beam  can 
be  obtained  within  a  reasonable  range  of  applied  current.  We  have  fabricated  and 
characterized  such  a  laser.  The  device  maintains  near  diffraction-limited  operation  up  to 
1.9  times  threshold  with  far  field  divergence  angle  less  than  0.3°  (FWHM).  A  surface- 
emitting  power  of  150  mW  is  obtained  under  continuous  wave  operation. 

Improvements  can  be  made  in  several  ways.  Using  bandgap-shifting  technique, 
we  can  design  the  grating  so  the  lasing  wavelength  is  at  the  gain  peak.  This  will  increase 
the  device  efficiency.  It  will  also  decrease  the  a-factor  and  suppress  filamentaion.  Better 
heatsink  and  packaging  can  be  used  to  decrease  the  thermal  effect  and  improve  the  device 
performance  at  higher  current. 
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Figure  Captions 


Figure  1 :  Schematic  top  view  of  a  grating  coupled  tapered  surface  emitting  laser. 

Figure  2:  Calculated  far  field  intensity  profiles  scattered  perpendicular  to  the  surface 

of  a  linear  grating. 

Figure  3:  Lateral  far  field  intensity  profile  for  a  circular  cuiwed  grating  at  1 .03  times 

the  threshold  current  Jn,. 

Figure  4:  Lateral  far  field  intensity  profile  for  a  circular  curved  grating  at  2.05  Jn,. 

Figure  5:  Lateral  far  field  intensity  profile  for  a  circular  curved  grating  at  3.08  Jn,. 

Figure  6:  Light-current  characteristics  of  the  surface-emitting  laser. 

Figure  7:  Far  field  of  the  surface-emitting  laser. 

(a)  I  =  1.1  A  (1.6  luO,  0FWHM  =  0.25°  (1.1  D.L.). 

(b)  1  =  1.3  A  (1.9  l,h),  0FWHM  =  0.30°  (1.4  D.L.). 
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